AJ, IN PRESS 

Preprint typeset using L^TgX style emulateapj v. 25/04/01 



M75, A GLOBULAR CLUSTER WITH A TRIMODAL HORIZONTAL BRANCH. 
I. COLOR-MAGNITUDE DIAGRAM 

M. Catelan 

Pontificia Universidad Catdlica de Chile, Departamento de Astronomfa y Astrofi'sica, Av. Vicuna Mackenna 4860, 690441 1 

Macul, Santiago, Chile 
mcatelan@astro.puc.cl 

J. BORISSOVA 

Institute of Astronomy, Bulgarian Academy of Sciences and Isaac Newton Institute of Chile Bulgarian Branch, 
72 Tsarigradsko chaussee, BG- 1784 Sofia, Bulgaria 
jura@ haemimont.bg 

F. R. FERRARO 

Osservatorio Astronomico di Bologna, Via Ranzani 1, 1-40127 Bologna, Italy 
ferraro @ apache . bo. astro it 

T. M. CORWIN 

Department of Physics, University of North Carolina at Charlotte, Charlotte, NC 28223 
mcorwin @ uncc.edu 

H. A. Smith 

Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824-1 116 

smith@pa.msu.edu 

AND 
R. KURTEV 

Department of Astronomy, Sofia University and Isaac Newton Institute of Chile Bulgarian Branch, James Bourchier Ave. 5, 

BG- 1164 Sofia, Bulgaria 
kurtev@phys.uni-sofia.bg 

AJ, in press 

ABSTRACT 

Deep U BVI photometry for a large field covering the distant globular cluster M75 (NGC 6864) is presented. We 
confirm a previous suggestion (Catelan et al. 1998a) that M75 possesses a bimodal horizontal branch (HB) bearing 
striking resemblance to the well-known case of NGC 1851. In addition, we detect a third, smaller grouping of 
stars on the M75 blue tail, separated from the bulk of the blue HB stars by a gap spanning about 0.5 mag in V. 
Such a group of stars may correspond to the upper part of a very extended, though thinly populated, blue tail. 
Thus M75 appears to have a trimodal HB. The presence of the "Grundahl jump" is verified using the broadband 
U filter. We explore the color-magnitude diagram of M75 with the purpose of deriving the cluster's fundamental 
parameters, and find a metallicity of [Fe/H] = -1 .03 ± 0.17 dex and -1 .24 ±0.21 in the Carre tta & Gratton (1997) 
and Zinn & West (1984) scales, respectively. We discuss earlier suggestions that the cluster has an anomalously 
low ratio of bright red giants to HB stars. A differential age analysis with respect to NGC 1851 suggests that the 
two clusters are essentially coeval. 

Subject headings: Stars: blue stragglers - Stars: Hertzsprung-Russell (HR) and C-M diagrams - Stars: 

horizontal-branch - Stars: variables: RR Lyr - Galaxy: globular clusters: individual: M75 
(NGC 6864) - Galaxy: globular clusters: individual: NGC 1851 



1. INTRODUCTION 

Adequate interpretation of the color-magnitude diagrams 
(CMDs) of globular clusters (GCs) is important for many as- 
trophysical and cosmological reasons. The horizontal branch 
(HB) evolutionary phase plays a particularly significant role in 
that regard. Because the observational properties of HB stars 
are extremely sensitive to the physical parameters of GCs, HB 
stars can impose invaluable constraints on the initial conditions 
that characterized the Galaxy in its infancy. 

Metallicity, [Fe/H], is the "first parameter" influencing HB 
morphology: Metal-rich clusters have, in general, more red HB 
stars than blue HB stars and vice versa for the metal-poor GCs. 
The existence of GCs which have the same [Fe/H] but different 
HB morphologies indicates that one or more second parameters 



exists. Despite the fact that the problem has received much at- 
tention since it was first discovered (Sandage & Wildey 1967), 
the second parameter(s) responsible for the spread in HB types 
at a given [Fe/H] has (have) not been firmly established. Under- 
standing the second parameter is arguably one of the most im- 
portant tasks that must be accomplished before the age and for- 
mation history of the Galaxy can be reliably determined (e.g., 
Searle & Zinn 1978; van den Bergh 1993; Zinn 1993; Lee, De- 
marque, & Zinn 1994; Majewski 1994). 

Much recent debate has focused on whether the predominant 
second parameter is age (e.g., Lee et al. 1994; Stetson, Van- 
denBerg, & Bolte 1996; Sarajedini, Chaboyer, & Demarque 
1997; Catelan 2000; Catelan, Ferraro, & Rood 2001b; Cate- 
lan et al. 2001a; VandenBerg 2000), but mass loss, helium 
and a-element abundances, rotation, deep mixing, binary in- 
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teractions, core concentration, and even planetary systems have 
been suggested as second parameters. It is likely that the second 
parameter phenomenon is in fact a combination of several in- 
dividual quantities and phenomena related to the formation and 
chemodynamical evolution of each individual star cluster, im- 
plying that the several individual second-parameter candidates 
have to be "weighted" on a case-by-case basis (e.g., Fusi Pecci 
et al. 1993; Fusi Pecci & Bellazzini 1997; Rich et al. 1997; 
Sweigart 1997a, 1997b; Kraft et al. 1998; Sweigart & Catelan 
1998; Soker 1998; Ferraro et al. 1999a; Catelan 2000; Soker & 
Harpaz 2000; Catelan et al. 2001a, 2001b). In this context, it 
has been extensively argued that an understanding of GCs pre- 
senting the second-parameter syndrome internally — i.e., those 
with bimodal HBs — is of paramount importance for under- 
standing the nature of the second-parameter phenomenon as a 
whole (e.g., Rood et al. 1993; Stetson et al. 1996; Borissova et 
al. 1997; Catelan et al. 1998a; Fusi Pecci & Bellazzini 1997; 
Walker 1998). 

Unfortunately, the current list of GCs known to have a bi- 
modal HB is not large; probably only NGC 1851, NGC 2808 
and NGC 6229 are universally recognized as "bona fide" 
bimodal-HB clusters, following the definition of the term pro- 
vided by Catelan et al. (1998a). However, as discussed by Cate- 
lan et al., there are several other GCs whose poorly investigated 
CMDs, once improved on the basis of modern CCD photome- 
try, may well move them into the "bimodal HB" category. Out- 
standing among those is M75 (NGC 6864). 

M75 is a moderately metal-rich ([Fe/H] ~ -1.3 in the Zinn 
& West 1984, hereafter ZW84, scale; Harris 1996), dense 
{log[/9 O /(M pc- 3 )] =4.9; Pryor & Meylan 1993} and dis- 
tant (Rq ~ 19 kpc) halo GC located about 13 kpc from the 
Galactic center (Harris 1996). Since it is not affected by 
high reddening — E(B-V) = 0.16 mag, according to the Harris 
catalogue — we find it somewhat surprising that the only CMD 
for this cluster available in the literature seems to be the pho- 
tographic one provided by Harris (1975). This is even more 
puzzling in view of the fact that it has been over a quarter of a 
century since Harris called attention to some intriguing charac- 
teristics of this cluster's CMD — such as an abnormally high ra- 
tio R of HB to red giant branch (RGB) stars and the presence of 
an unusually strong blue HB component for such a moderately 
metal-rich cluster. Interest in the CMD of this particular cluster 
seems to have been revived only very recently, when Catelan et 
al. (1998a) pointed out that "new CCD investigations are es- 
pecially encouraged [for M75], since the latest CMD available 
(Harris 1975) is strongly suggestive of an NGC 1851 -like HB 
bimodality." Note that, due to M75's distance, the Harris pho- 
tographic CMD was only barely able to reach the HB level of 
the cluster. 

It is the main purpose of the present paper to present the first 
deep CCD photometry of M75, and to investigate whether — as 
suggested by Catelan et al. (1998a) — M75 can be classified as 
yet another bona fide bimodal-HB cluster. In §2, we discuss our 
observing runs and data reduction techniques. In §3, the cluster 
CMD is discussed in detail, revealing not only a clearly bimodal 
HB but also a third grouping of blue tail stars clearly separated 
from the bulk of the blue HB, and likely indicating the cooler 
end of a very extended, though thinly populated, blue tail. The 
presence and radial distribution of blue straggler stars (BSS) is 
also discussed, and Harris's (1975) suggestion of an abnormally 
high R ratio is addressed. In §4, we evaluate the "photometric 
metallicity" and reddening of the cluster, and in §5 we derive 



its age relative to NGC 1851. 

In a companion paper (Corwin et al. 2001, hereinafter Pa- 
per II), we will present a detailed analysis of the variable star 
population in the cluster. 

2. OBSERVATIONS AND DATA REDUCTION 

Our analysis is based on the following observational mate- 
rial: 

1. A set of "blue" U,B and "red" V,I CCD frames ob- 
tained in June 1997 on the ESO-NTT telescope with the 1024 x 
1024 and 2048 x 1900 CCD cameras. The array scales were 
OJe^'pixel" 1 and 0.266" pixel" 1 , giving fields of view of about 
6.16' x 6.16' and 9.08' x 8.42'. The data consist of one short 
(60-20 sec) exposure image in B,V,I, and two long exposure 
images (900-300 sec) in U,B,V,I. 

2. Fourteen B,V frames obtained in July 1999 at the 0.9- 
meter telescope of the Cerro Tololo Inter-American Observa- 
tory with a 2048 x 2048 CCD camera. The pixel scale of 
0.396"pixel _1 gave an observing area of 13.5' x 13.5'. 

3. Ten B,V frames obtained in July 1989 at the 1.54-meter 
Danish telescope with a 512 x 320 CCD camera. The pixel 
scale of 0.47" pixel -1 gave an observing area of 4' x 2.5'. 

The stellar photometry of NTT images was carried out sepa- 
rately for all frames using daophot/allstar (Stetson 1993). 
The instrumental magnitudes are transformed photometrically 
to the deepest reference frame in each filter. The magnitudes 
in the same bands of the common unsaturated stars are av- 
eraged. The brightest stars are measured only on the short 
exposure time images. The field of view of "blue" (U and 
B) images is smaller than for "red" (V and /) images. The 
datasets were matched together and a final catalog with stel- 
lar coordinates and the instrumental magnitudes in each filter 
has been compiled for all objects identified in each field. The 
instrumental values of all the NTT frames were transformed 
to the standard system using 20 standard stars in five selected 
fields from Landolt (1992), covering a color range -0.24 < 
(B-V) < 1.13. The mean photometric errors are w 0.04 mag 
for U, B,V,I < 22 mag and between 0.06 mag and 0.1 mag 
for the fainter magnitudes. The formal errors for all stars vs. 
their magnitudes are displayed in Fig. 1. The artificial star 
technique (Stetson & Harris 1988; Stetson 1991a, 1991b) was 
used to determine the completeness limits of the data. The 
50% completeness limit of the photometry is determined at 
(U,B,V,I) = (22, 22.4, 22.4, 21 .8) mag at a distance of 0.8' from 
the cluster center. 

daophot/allframe (Stetson 1993) was applied to 28 of 
the CTIO images (14 B and 14 V). All frames are calibrated 
independently using 1 3 photoelectric standards from Alvarado 
et al. (1990) and Harris (1975). The photoelectric standards 
included stars as red as the reddest M75 red giant stars. A file 
with the coordinates of 7504 stars and 28 columns of standard 
magnitudes (14 columns of V and 14 of B-V) was obtained. 
For each star the mean V magnitude and B-V color as well as 
the standard deviations were determined. The mean photomet- 
ric errors are w 0.08 mag for B, V < 19 mag and 0.12 mag for 
the fainter magnitudes. The 50% limit of the photometry is at 
(B,V) = (21,21.2) mag at radius 1'. 

The stellar photometry of 1.54m Danish images was also 
carried out with DAOPHOT/ALLSTAR. The mean photometric 
errors are « 0.03 mag for B, V < 19 mag and 0.08 mag for 
the fainter magnitudes. The 50% limit of the photometry is at 
(B,V) = (18.5, 18.3) mag at radius 0.3'. The seeing for the Dan- 




FIG. 1 . — Formal errors of the NTT photometry from the DAOPHOT package as a function of the mean U, B. V. I magnitudes for each star. 



ish observing run was better (0.9") than for the NTT observing 
run (1.2"). 

We transformed the CTIO and Danish x,y coordinate systems 
to the NTT local system and then searched for stars in common 
between the two datasets. We checked the residuals in mag- 
nitude and color for the stars in common. They do not show 
any systematic difference or trend so we can conclude that all 
datasets are homogeneous in magnitude within the errors. 

Taking into account the above errors and completeness de- 
termination, we used the NTT dataset for the global analysis 
of the CMD. For star counts we added 28 bright stars from the 
Danish dataset which are saturated on NTT images and have 
radial distances between 0.3' and 0.8'. It should be noted that 
these 28 stars were used only for number counts. The fiducial 
line of M75 (see §3.1 below) was determined using solely the 
NTT dataset. 

The field stars are statistically decontaminated using off-field 
CTIO images taken 30' east from the center of the cluster and 
normalized to the same area as the NTT field of view. The 
(B-V,V) CMDs of "cluster + field" and "field" are divided into 
15 boxes (shown in Fig. 2) — five intervals in magnitude and 
three in color. The stars in each box in the two diagrams are 



counted. Then, an equivalent number of stars is removed from 
the single boxes of the "cluster + field" CMD on the basis of 
the number of field stars found in the "field" CMD alone (right- 
hand-side panel in Fig. 2). 

3. THE COLOR-MAGNITUDE DIAGRAM 
3.1. Overall CMD Morphology 

The (V, B-V) and (V, B-I) CMDs from the NTT dataset are 
presented in Fig. 3. Field stars were statistically subtracted as 
described in §2, and known variables are omitted from the plot. 

Mean ridgelines of the main branches of M75 were deter- 
mined from the (V, V-I) CMD. With this purpose, the fol- 
lowing selection criteria were adopted in order to have the best 
definition of the CMD branches: for V < 19.7 mag, all stars 
were used; for V > 19.7 mag, and in order to better outline 
the turnoff (TO) region, only stars with r > 150" were em- 
ployed (thus avoiding the crowding in the innermost regions 
while minimizing field star contamination). 

The mean ridgelines of the main branches in the CMD were 
determined following the polynomial fitting technique by Sara- 
jedini (1994). A first rough selection of the candidate RGB 
stars was performed by eye, removing the HB and part of the 
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FIG. 2. — (B-V.V) color-magnitude diagram for M75 (left panel). The field CMD (right panel) is based on images taken 30' east of the center of the cluster. 



asymptotic giant branch (AGB) stars. A polynomial law in the 
form (V —I) = F(V) was used. In several iterations the stars 
deviating by more than 2-sigma in color from the fitting rela- 
tion were rejected. The standard error of the resulting fit is 
0.04 mag. The mean ridgelines of the fainter part of the CMD 
(V > 19.7 mag — subgiant branch and main sequence) were de- 
termined by dividing these branches into bins and computing in 
each bin the mode of the distribution in color. The "ridgeline" 
for the red HB region was determined as the lower envelope of 
the red clump. Table 1 presents the adopted normal points for 
each branch. In Fig. 4 the (V, V-I) CMD is shown with the 
mean ridgelines overplotted. 

The main sequence (MS) TO point is found to be at V T o = 
21.22 ±0.09 mag and (V-I) T0 = 0.818 ±0.05 mag. The mag- 
nitude level of the HB, determined as the lower boundary of 
the red HB "clump," is V HB = 17.80 ± 0.03 mag. This im- 
plies a magnitude difference between the HB and the TO of 
AV™ = 3. 42 ±0.09 mag. 

3.2. Metallicity and Reddening 

The best fit of (B-V) vs. (U-B), (B-V) vs. (V-I), and 
(U-B) vs. (V-I) distributions for the M75 blue HB stars 
to the color-color lines from Bessel (1990) gives E(B-V) = 
0. 16 ± 0.02. The fit is shown on Fig. 5. 

For comparison, the Schlegel, Finkbeiner, & Davis (1998) 
COBE/DIRBE dust maps give E(B-V) = 0.152 mag at the M75 
position, and Harris (1996) lists E(B-V) = 0.16 mag in the lat- 
est (June 1999) version of his catalogue. 



A complete set of RGB parameters and metallicity indica- 
tors in the classical (V, B-V) plane has been recently presented 
by Ferraro et al. (1999a, hereafter F99). F99 also reported 
on an independent calibration of RGB parameters in terms of 
the cluster metallicity (both in the Carretta & Gratton 1997, 
hereafter CG97, [Fe/H] scale and in the corresponding global 
[M/H] scale). In particular, F99 presented a system of equa- 
tions (see Table 4 in F99) which can be used to simultaneously 
derive an estimate of metal abundance (in terms of [Fe/H] CG97 
and [M/H]) and reddening from the morphology and location 
of the RGB in the (V, B-V) CMD. The "global metallicity" 
[M/H] assumes an a-element enhancement of [a/Fe] = 0.28 
for [Fe/H] <-l. 

For M75, we used the mean ridgeline to measure the follow- 
ing RGB parameters: (B-V)o tg , defined as the RGB color at 
the HB level; the two RGB slopes S2.0 an d S2.5, defined as the 
slope of the line connecting the intersection of the RGB and 
HB with the points along the RGB located, respectively, 2.0 
and 2.5 mag brighter than the HB; AVi.i, AV L2 and AV1.4— 
defined as magnitude differences between the HB and RGB 
at the fixed colors (B-V)o = 1.1,1.2 and 1.4 mag. The de- 
rived metallicity, based on the weighted mean of the above de- 
fined values as listed in Table 2, is [Fe/H] CG97 = -1.01 ±0.12, 
[M/H] = -0.82 ±0.19 and [Fe/H] zw84 = -1.22 ±0.21. The 
derived reddening is E(B-V) = 0.155 ± 0.032. The errors as- 
sociated with the measurements are conservative estimates of 
the global uncertainties, formal errors being much smaller than 
those assumed. 
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FIG. 3.— (V, B-V) and (V, B-I) CMDs for all stars in the NTT dataset for M75. Known variable stars are omitted from the plot. 



3.3. Red Giant Branch 

The so-called RGB "bump," corresponding to the point in 
the evolution of RGB stars when the H-burning shell encoun- 
ters the chemical composition discontinuity left behind by the 
maximum inward penetration of the convective envelope, was 
first predicted by theoretical models (Thomas 1967; Iben 1968) 
and subsequently observed in many clusters (see F99 for more 
details). The break in the slope of the cumulative RGB lumi- 
nosity function is a common technique used to locate the RGB 
bump (Fusi Pecci et al. 1990). The observed differential and 
integrated luminosity functions for the M75 RGB are shown in 
Fig. 6. As can be seen the RGB bump of M75 is very well 
defined at V = 17.75 ± 0.05 mag. 

The magnitude difference between the RGB luminosity func- 
tion "bump" and the HB level, AV^g mp , is a metallicity indicator 
recently recalibrated by F99. Using the data given in their Ta- 
ble 5 and the calibration equations from Table 6 in F99 we cal- 



culated [Fe/H] CG97 =-1.10 ±0.10, [M/H] =-0.91 ±0.09 and 
[Fe/H] zw84 = -1.30 ±0.08. As can be seen, the RGB bump- 
based metallicity is in good agreement with the value derived 
in §3.2. A plot showing the difference in brightness between the 
RGB bump and the HB, AV^ P , as a function of metallicity is 
presented in Fig. 7. 

Taking the RGB bump-based metallicity into account, the 
revised metallicity of M75, based on the weighted mean over 
all the values listed in Table 2, is [Fe/H] CG97 = -1.03 ±0.17 
and [M/H] = -0.87 ±0.19 on the Carretta & Gratton (1997) 
scale. On the Zinn & West (1984) scale we find [Fe/H] zw84 = 

1.24 : 0.21. 

Recently Bono et al. (2001) defined the new parameter 
^bump as the ratio between the number of stars in the bump 
region (Vb um p ± 0.4) and the number of RGB stars in the in- 
terval Vbump + 0.5 < V < Vbump + 1-5. In our case we calculated 
flbum P = 0.47 ±0.04. 
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FIG. 4. — (V, V -I) CMD for M75, with the derived ridgeline overplotted. Crosses indicate variable stars (from Paper II). The dashed line indicates the magnitude 
level which separates the two samples used for better definition of the CMD branches and determination of the ridgeline: while for V < 19.7 mag all stars were 
utilized (and are displayed), fainter than this level only stars with r > 150" were employed (and are shown). Note that, for the red HB, the "ridgeline" actually 
represents the lower envelope of the distribution. 



Comparison of the Rb um p parameter with the cluster sample 
of Bono et al. (2001) — their Table 1 — shows good agreement. 
In Fig. 8, where /?bump is plotted as a function of the metallicity, 
M75 is marked with a filled circle, showing its normal location 
with respect to the other clusters. This "normal" behavior of 
^bump in M75, in analogy with the discussion presented by Bono 
et al. for the other clusters plotted in Fig. 8, suggests that the 
M75 bump too has not been affected by the occurrence of non- 
canonical deep mixing that might have substantially altered the 
chemical composition profile in the vicinity of the H-burning 
shell in the early evolution of M75's red giants. 

Harris (1975) first mentioned an abnormally high ratio R of 
HB to red giant branch (RGB) stars. To calculate the R ratio we 
find the following number counts: A^hb = 149, Nrgb = 118, and 
A^agb = 28 — where Nrgb is the number of RGB stars brighter 
than the HB level, following the definition by Buzzoni et al. 
(1983). The observed star counts are corrected for complete- 
ness and the estimated numbers rounded to the closest inte- 
ger value. It thus follows that R = N H b/N RG b = 1.26 ±0.15. 
For R' = AWN(rgb+agb), the value is R' = 1.02 ± 0.10. For 
Rl = N agb /N RC b we calculated Rl = 0.24 ± 0.05 and R2 = 
Nagb/Nhb is Rl = 0.19 ± 0.04. We have assumed, following 
Buzzoni et al., a differential bolometric correction between the 
HB and the RGB of 0. 15 mag. 



If we adopt the Zoccali et al. (2001) definition of jVrgb as 
RGB stars brighter than the ZAHB (which is approximately 
equivalent to the lower envelope of the red HB), A^rgb be- 
comes 90 (after completeness correction) and R = 1 .66 ± 0.20, 
R' = 1 .26 ± 0. 16 andtfl = 0.3 1 ± 0.07. Plotting the R parameter 
vs. global metallicity in Fig. 9 (the cluster data are from Table 1 
in Zoccali et al.) we can see that M75 has a high R value for its 
metallicity — though only at the 1 - 2a level. 

Note that the R value following the Buzzoni et al. (1983) 
definition does not seem atypical (or large). It is only when we 
use the Zoccali et al. (2001) definition that we are able to verify 
the suggestion by Harris (1975) of a low ratio of giants to HB 
stars in M75. This problem is due to the fact that the bump, at 
such a metallicity, is critically located close to the HB level (see 
Fig. 7), and may or may not be included in the number counts 
due to only a small change in the magnitude cutoff. 

In order to further investigate this issue, we have utilized 
the observed RGB number counts down to a deeper magni- 
tude level on the CMD, performing a direct comparison against 
similarly computed values from the theoretical simulations dis- 
cussed in §3.4.1 and presented in Figs. 11 and 12 below. In 
order to avoid the ambiguity brought about by the position- 
ing of the RGB bump (see Fig. 7) in the derived number 
counts, we define i?(M„)+i .o an d R{m v )+i.5 as the ratio between 




FIG. 5. — Dereddened color-color (B-V)q vs. (U-B)q, (B-V)q vs. (V-I)q, and (U-B)o vs. (V-/)o diagrams. The dashed lines are sequences of luminosity class 
V from Bessel (1990). The arrows show the direction of the reddening vector. 



HB stars and RGB stars brighter than (M y (HB)} + 1 mag and 
(My(HB)) + 1.5 mag, respectively. The resulting comparison 
between the observed and theoretical number ratios is given in 
Table 3. As can be seen, the observed R -ratios are clearly con- 
firmed to be higher than the values predicted for T M s = 0.23, 
in agreement with the conclusion based on the Zoccali et al. 
(2001) definition of the R -ratio and with the original sugges- 
tion by Harris (1975) that M75 seems to have too few giants 
compared to HB stars. 

What is the origin of the larger 7?-ratio for M75? The tradi- 
tional explanation would involve a larger helium abundance for 
this cluster (Iben 1968). Is this a viable explanation, in the case 
of M75? In order to answer this question, we may take bene- 
fit of the predicted dependence of HB luminosity on the helium 
abundance, which should accordingly lead to longer periods for 
M75's RR Lyrae variables (e.g., Catelan 1996). We will check 
whether this is the case, or whether other explanations are re- 
quired, in Paper II, where we present a detailed analysis of the 
M75 RR Lyrae light curves. 

3.4. HB Morphology 

3 .4. 1 . HB Multimodality 

Fig. 10 shows the HB region of the CMD of M75 in sev- 
eral different planes: (V, U -V); (V, U-I); and (U, U-V). 
Known variable stars (see Paper II) are omitted from these 
plots. Both the red and the blue HB regions are very well pop- 
ulated, as first noted by Harris (1975). The blue part of the 
HB has a blue tail and an obvious "gap" at V ~ 18.5 mag and 
(U-V)o between -0.55 and -0.15, corresponding to tempera- 



tures in the range 10,500 < T e ff(K) < 13,500 — which overlaps 
the suggested "Gl" and "G2" gaps in Ferraro et al. (1998). An- 
other, much less obvious gap (or underpopulated region) may 
be present at V « 17.9- 18.0 mag and (U-V) w 0.1 mag, 
(B-/)o«0.1mag. 

In order to more quantitatively describe the morphology of 
the HB, we computed several useful HB morphology parame- 
ters, including those defined by Mironov (1972), Zinn (1986), 
Lee, Demarque, & Zinn (1990), Buonanno (1993), Fusi Pecci 
et al. (1993), and Catelan et al. (2001b). Several of these au- 
thors have emphasized the importance of utilizing as many HB 
morphology parameters as possible when drawing general in- 
ferences based on the HB morphology of Galactic GCs. The 
values of the computed parameters are given in Table 4. We 
briefly recall that B, V, and R are the number of blue HB, vari- 
able (RR Lyrae), and red HB stars, respectively; B2 is the num- 
ber of blue HB stars bluer than (B-V) Q = -0.02 mag; B2' is the 
number of HB stars redder than (V— 7)o = -0.02 mag; BT is 
the number of blue HB stars redder than (V-I)o = +0.07 mag; 
and, finally, B0' is the number of blue HB stars redder than 
(y-/) = 0.0 mag. 

To compute the reddening-dependent HB morphology pa- 
rameters, we have assumed E(B-V) = 0.16 (§3). Also, the 
number of variable stars, V = 22, was adopted from our compre- 
hensive variability survey, described in Paper II. All observed 
number counts are corrected for completeness. 

It is obvious from these numbers that M75 is indeed another 
cluster with a strongly bimodal HB, since it has many fewer 
RR Lyrae variables than either red HB or blue HB stars. A 
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FIG. 6. — The observed differential (upper panel) and integrated (lower panel) luminosity function for the M75 RGB. 
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histogram with the (U - V)o color distribution of the M75 HB 
stars (corrected for completeness) is given in the bottom panel 
of Fig. 10. Since we do not have U — V mean colors for the 
variables we assume a flat distribution. This clearly shows that 
the distribution of stars along the HB of M75 has three main 
modes — though the third mode, on the blue tail, is very thinly 
populated. For a detailed discussion of the completeness of our 
RR Lyrae search, see Paper II. 

The M75 HB clearly bears striking resemblance to that of 
the well-known bimodal-HB GC NGC 1851, as first pointed 
out by Catelan et al. (1998a) on the basis of the Harris (1975) 
photographic CMD. In particular, the number ratios B : V : R = 
0.38 : 0.16 : 0.45 are similar to those for NGC 1851, for which 
B : V : R = 030 : 0.10: 0.60. 

We have investigated further the issue of HB bimodality in 
M75 by computing a new grid of synthetic HBs aimed at match- 
ing the HB morphology parameters displayed in Table 4. Our 
purpose is to check whether a bimodality in the underlying 
evolutionary parameters is necessarily implied by the observed 
bimodality in the HB number counts (Catelan et al. 1998a). 
Our simulations, computed with SINTDELPHI (Catelan et al. 
2001b), utilize the evolutionary tracks for Y MS = 0.23, Z = 0.002 
from Sweigart & Catelan (1998). Observational scatter has 



been added by means of suitable analytical representations of 
the errors in the photometry. The number of HB stars in these 
simulations is similar to the observed one, Nub = 149, but has 
been allowed to fluctuate according to the Poisson distribution. 

We first attempted to reproduce the observed HB morphol- 
ogy parameters for M75 by means of a unimodal Gaussian de- 
viate in the zero-age HB (ZAHB) mass. Following a procedure 
similar to that described in Catelan et al. (2001a, 2001b), we 
find a best match between the observed and predicted HB mor- 
phology parameters for (M H b) = O.624M , a M = 0.043M Q . 
Some randomly selected samples of simulations for this uni- 
modal mass distribution case are displayed, in the [M v , (B-I)o] 
plane, in Fig. 11. The predicted HB morphology parameters for 
this case are shown in Table 4 (third column); the indicated "er- 
ror bars" represent the standard deviation of the mean over a set 
of 1000 simulations. It is readily apparent that, while the overall 
general agreement, particularly in terms of the commonly em- 
ployed Mironov and Lee-Zinn parameters, is reasonable, some 
of the other parameters introduced by Catelan et al., especially 
BO'/B and BT /(B + V + R), are not accounted for in an entirely 
satisfactory way by these simulations. For this reason, we have 
attempted to compute extra sets of simulations for an intrinsi- 
cally bimodal distribution in mass on the ZAHB. 
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In the bimodal case, we have found that the following com- 
bination of parameters gives the best fit between the models 
and the observations: (M H b i) = 0.595 M Q , a M \ = 0.010M Q , 
JVhb,i = 46; (M HB ,2) = 0.640 M Q , a M , 2 = 0.030M Q , N w ,2 = 104. 
Some randomly selected samples of simulations for this bi- 
modal mass distribution case are displayed, in the [M v , (B-I)o] 
plane, in Fig. 12. The predicted HB morphology parameters for 
this case are also shown in Table 4 (fourth column). 

We find that the intrinsically bimodal mass distribution case 
provides a much more satisfactory fit to the whole set of HB 
morphology parameters than does the unimodal case. To some 
extent, this is not an unexpected result, given that this case con- 
tains five free parameters, as opposed to two in the unimodal 
one (Catelan et al. 1998a). The main difference between the 
bimodal and unimodal cases, as can be readily appreciated by 
inspection of Figs. 1 1 and 12, is that the unimodal case clearly 
predicts that a fairly extended blue tail should be present in 
M75, whereas this is not necessarily required in the bimodal 
case. Also in this respect, apart from a few cases where sta- 
tistical fluctuations change the predicted simulations quite sub- 
stantially, the bimodal case appears much more successful at 



matching the observed distribution, showing, in general, a more 
stubby blue HB with relatively few hot stars being in general 
expected. Note, however, that neither the unimodal nor the bi- 
modal cases shown in Figs. 1 1 and 12 can successfully account 
for the existence of a third mode of hotter HB stars along the 
blue tail of M75 (see also Fig. 10, bottom panel). 

3.4.2. The Grundahl Jump in Broadband U 

Grundahl et al. (1999) have recently shown that the inter- 
pretation of Stromgren u photometry of blue HB stars in GCs 
is affected by a most intriguing phenomenon: at temperatures 
higher than T e ff « 1 1 ,500 K, all GC CMDs show an intriguing 
"jump" in the u, (u— y)o plane, with all blue HB stars hotter than 
1 1,500 K appearing brighter and/or hotter than predicted by the 
models. As argued by Grundahl et al., and spectacularly con- 
firmed by Behr et al. (1999), the reason for this discrepancy can 
be attributed to the fact that the atmospheres of hot HB stars in 
GCs are affected by radiative levitation of metals, which pumps 
up heavy metals to supersolar levels in these stars' atmospheres. 

As can be seen in Fig. 10 (bottom panel), a similar jump is 
present also in broadband U, with all stars blueward than (U - 
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V)o = -0.5 falling systematically above the canonical ZAHB. 
The displayed ZAHB sequence corresponds to VandenBerg et 
al. (2000) models for [Fe/H] = -1.14 and [a/Fe] = +0.3, and 
were transformed to the observational plane using standard pre- 
scriptions from Kurucz (1992) — which assume that the clus- 
ter metallicity is appropriate also for the atmospheres of hot- 
ter blue HB stars. It is worth noting that the Grundahl jump 
phenomenon in broadband U has recently been verified also by 
Siegel et al. (1999), Bedin et al. (2000), and Markov, Spassova, 
& Baev (2001). The presence of the jump phenomenon in 



broadband U is in good agreement with the radiative levitation 
scenario, as can be seen from inspection of Fig. 7 in Grundahl 
etal. (1999). 

3.5. Blue Straggler Stars 

On the basis of their location in the (V, U-I) CMD, we iden- 
tified a sample of 26 candidate BSS. The x and y coordinates (in 
pixels) and U, B, V, I magnitudes and distance from the cluster 
center (in arcmin) for these candidate BSS are listed in Table 5. 

Fig. 13 shows zoomed M75 CMDs where the BSS candidates 
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are plotted as encircled dots. As can be seen, all the candidates 
selected in the (V, U-I) plane lie also in the BSS region in 
the (U, U-I) CMD. Moreover, the latter CMD suggests that a 
few additional objects might be considered BSS candidates. We 
conservatively decided, however, to limit our analysis to those 
stars which are BSS candidates in both the indicated CMDs. 

Careful checks were made to ensure that the internal errors of 
the BSS candidates are the same as those of the subgiant branch 
(SGB) stars at the same brightness level. The DAOPHOT param- 
eters SHARP and CHI were checked for each BSS candidate, in 
order to ensure that most are unlikely to be the result of spuri- 
ous photometric blends (e.g., Ferraro, Fusi Pecci, & Buonanno 
1992). The positions of the BSS were checked on the statisti- 
cally decontaminated CMDs. The number of field stars is not 
sufficient to explain the BSS population (see also Fig. 2). 

In order to quantitatively compare the BSS populations 
among different Galactic GCs, one needs to suitably normalize 



the numbers of such stars. In this sense, Ferraro et al. (1999b) 
defined the specific fraction of BSS as the number of BSS nor- 
malized to the number of HB stars observed in the same cluster 
region, F^ s = Nbss/Nhb- In M75, the specific frequency cal- 
culated in this way is F§^ s ~ 0.15. 

Generally, BSS are more centrally concentrated with respect 
to the other stars in the cluster, although there are notable ex- 
ceptions to this general rule, such as M3 (Ferraro et al. 1993, 
1997) and M13 (Paltrinieri et al. 1998). In M75 we compared 
the radial distribution of BSS with respect to SGB stars (in the 
range 19.5 < V < 20.5). The cumulative radial distributions 
are plotted, as a function of the projected distance (r) from the 
cluster center, in Fig. 14, where one sees that the M75 BSS 
(circles) are more centrally concentrated than the SGB sam- 
ple (diamonds). The Kolmogorov-Smirnov test shows that the 
probability of drawing the two populations from the same dis- 
tribution is 0.3%. 
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FIG. 1 1. — Random sample of HB simulations for M75 in the unimodal ZAHB mass distribution case. The My, (B—])g plane is shown. We have assumed 
E(B—1)/E(B-V) = 2.7, after Stetson (1998). In the upper left panel, the observational data are shown, along with the reddening value and distance modulus found 
more suitable to fit the theoretical ZAHB for the indicated chemical composition to the observed distribution. Variable stars are omitted in the panel showing the 
observed CMD, but not in the others (where they are displayed as encircled gray dots) — which show the theoretical CMD simulations. Also indicated are the 
B : V : R number ratios and the value of the Buonanno (1993) HB morphology parameter. 



4. AGE DETERMINATION 

Two main methods are currently used to measure relative 
ages: (i) the so-called vertical method — based on the TO lumi- 
nosity with respect to the HB level (Buonanno, Corsi, & Fusi 
Pecci 1989) and (ii) the horizontal method — based on the accu- 
rate determination of the color difference A(B-V)jq B between 
the base of the RGB and the TO (VandenBerg, Bolte, & Stetson 
1990). 

Admittedly, our photometry in the [(B—V), V] plane is not 
sufficiently accurate to derive precise color difference measures 
down to the TO region. In addition, the horizontal method is 
so intrinsically sensitive to small variations in colour that even 
small (<~ 0.01 mag) uncertainties in the registration of the mean 



ridge lines could generate large (~ 1 Gyr), "artificial" differ- 
ences in the derived ages. For this reason, in order to obtain a 
relative age for M75, we applied only the vertical method. 

In doing this, we perform a comparison with respect to an- 
other cluster with similar metallicity and HB morphology — 
namely, NGC 1851. Recent photometry in the [(V-I), V] plane 
has been published for NGC 1851 by Bellazzini et al. (2001). In 
Fig. 15 the fiducial lines of NGC 1851 (from Bellazzini et al.) 
are overplotted on the CMD of M75 in the [(V-I) , V] plane. 

As can be seen from Fig. 15, while both the HB and the 
SGB/TO regions of the two clusters nicely agree, a small dif- 
ference in colour in the brightest region of the RGB does exist 
between the two data sets. This small colour residual could in 
principle be due to a small difference in the global metallicity 



13 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



Ml I | I I I I | I I I 



I 1 1 1 1 I 1 

» V K'3M DM □*? 








-1 

o 
1 

2 
3 



■ UmtuUlhlfelllMt'Y'Oqj'DDIi; 
- » V K'DSr Dib D4S f* 



1 

2 



> ^ 



/ :: 




"» 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 i ' 

4i m IS 1..J fc* ±4 2.S 4i 44 &3 1.0 1.4 2j* 



■- WittuUihiakii 

V ft' El 3b DM □ 



m.Y-DHI-DDDJ 




» V R'DS DID c« 

P-7 FYP-'V'*:i ■" ■□!] 




FIG. 12. — As in Fig. 1 1 , but for the bimodal mass distribution case. 



between the two clusters, although the metallicity value derived 
in §3.2 for M75 is very similar to that from Ferraro et al. (1999) 
for NGC 1851. In our opinion, the noted mismatch is likely due 
to a not completely accounted for "colour equation" in one of 
the two datasets. This effect is generally due to poor coverage 
in colour of the standard stars. 

On the basis of the VandenBerg & Bell (1985) models, Buo- 
nanno et al. (1993) obtained the following "vertical" relation 
between the differential age parameter A = Ai V™ - A2V™ 
and the age tg in Gyr: 



Alogf 9 = (0.44 + 0.04[Fe/H])A. (1) 

In our case, we shall assume that subscript "1" stands for M75, 
and "2" for NGC 1851. For NGC 1851 we used the data given 
in Bellazzini et al. (2001) to find AV T H B = 3.43 ±0.1. Assum- 
ing [Fe/H] ~ -1.05 for both M75 (§3.2) and NGC 1851 and 
using the value of AV^ derived in §3.1 for M75, we obtain 
A = 3.42-3.43 = -0.01 ±0.1 mag, implying Alogf 9 - -0.004 
and thus 
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A f 9 (M75 - NGC 1851) 0.1 ±2.0 Gyr. 

This result suggests that M75 is essentially coeval with 
NGC 1851. For recent comparisons of the age of NGC 1851 
with that of other GCs with similar metallicity, including 
NGC 288 and NGC 362, the reader is referred to Rosenberg 
et al. (1999), VandenBerg (2000), and Bellazzini et al. (2001); 
these papers all seem to indicate that NGC 1851 is (slightly) 
younger than other GCs of similar metallicity. 

We note that even though the Buonanno et al. (1993) for- 
mula, Eq. ( 1 ) above, is not based on the latest evolutionary mod- 
els available in the literature, they are still expected to give suf- 
ficiently accurate relative ages, under the hypothesis that both 
clusters have the same chemical composition. 

5. SUMMARY AND CONCLUDING REMARKS 

In this paper, we have provided the first CCD photometry for 
the distant GC M75. We confirm a previous suggestion (Cate- 
lan et al. 1998a) that this cluster contains two main modes on 
the HB, and show that the cluster actually has a third mode on 
the blue tail, clearly separated from the bulk of the blue HB 



by a wide gap, encompassing about 3000 K in temperature. In 
CMDs utilizing broadand U magnitudes, these hotter HB stars 
are found to be clearly inconsistent with canonical evolution- 
ary models, showing that the "Grundahl jump" phenomenon 
(Grundahl et al. 1999) can be studied using broadband U be- 
sides Strtimgren u. We also find a population of centrally con- 
centrated BSS in M75. The metallicity inferred for this cluster 
from our photometry, [Fe/H] = -1.03 ±0.17 in the CG97 scale 
and [Fe/H] = -1 .24 ± 0.21 in the ZW84 scale, is similar to that 
of the well-known bimodal-HB globular NGC 1851, thus en- 
abling a straightforward comparison of their turnoff ages. We 
find that the two clusters have essentially the same age. This is 
an interesting result, in view of the fact that NGC 1851 itself 
has been argued to be slightly younger than other GCs of simi- 
lar metallicity (Bellazzini et al. 2001 and references therein). 
However, better photometry is needed to conclusive demon- 
strate that M75 too is younger than the average for its metal- 
licity. 

What is the physical reason for the HB multimodality in 
M75? At present, this remains unclear. Similar anomalies in 
HB morphology have been seen in other clusters, and the sug- 
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gested explanations — none of which is seemingly capable of 
accounting for all the peculiarities — usually involve environ- 
mental effects, stellar rotation, abundance anomalies, mass loss 
in RGB stars, binarity, and even planets (e.g., Sosin et al. 1997 
and references therein). In the case of M75, we are in a bad po- 
sition to analyze most of these factors, given the lack of spec- 
troscopic data for the cluster stars. 

However, at least as far as the HB bimodality (i.e., relative 
absence of RR Lyrae stars, compared with both blue and red 
HB stars) is concerned, we can make use of RR Lyrae pulsa- 
tion properties in order to investigate the presence of any non- 
canonical effects that might impact their pulsational properties 
in a significant way. This could at least differentiate processes 
which affect HB morphology by means solely of RGB mass 
loss from those that may affect not only the HB star masses, but 
also their luminosities (Catelan, Sweigart, & Borissova 1998b). 
Similarly, useful constraints can thereby be placed on the origin 
of M75's peculiar R -ratio. 

With this in mind, we undertook an extensive variability sur- 



vey of M75. In Paper II we will provide a detailed discussion of 
the variable stars population in this cluster, in the hope that this 
will help us better understand the possible origin(s) of its pecu- 
liar HB morphology. Accordingly, a more extensive discussion 
of possible explanations for M75's peculiar HB will also be de- 
ferred to Paper II. 
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Table 2 
M75 Metallicity 



Parameter [Fe/H] 



CG97 scale 




(fi-V)og -1.21 ±0.22 
AVu -0.95 ±0.16 
AVi 2 -0.92 ±0.14 
AV14 -0.97 ±0.12 
S 25 -0.74 ±0.19 
S 2 .q -0.74 ±0.19 
AV^ mp -1.10 ±0.10 




ZW84 scale 




(V-B)o.g -1.29±0.15 
AVu -1.21 ±0.14 
AVi.2 -1.16 ±0.21 
Ay HB mP -1.30±0.13 




Table 3 

Ratio between HB and RGB Stars in M75 


Parameter Observed Model, Unimodal 


Model, Bimodal 


R {Mv ) 1.75 ±0.20 1.52 ±0.20 
R {Mv )+\o 0.89±0.12 0.73 ±0.08 
R(m v ) + \5 0.80±0.09 0.51 ±0.05 


1.52 ±0.20 
0.73 ±0.08 
0.51 ±0.05 
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Table 4 

HB Morphology Parameters for M75 



Parameter Observed Model, Unimodal Model, Bimodal 



B:V:R 


0.38: 0.16: 0.45 


0.369: 0.164: 0.467 
(0.040) (0.030) (0.040) 


0.374: 0.162: 0.463 
(0.040) (0.030) (0.039) 


B/{B + R) 


+0.46 ±0.05 


+0.442 ±0.045 


+0.447 ±0.044 


(B-R)/(B + V+R) 


-0.07 ±0.03 


-0.097 ±0.075 


-0.089 ±0.073 


(B2-R)/(B + V+R) 


-0.29 ±0.05 


-0.251 ±0.064 


-0.332 ±0.056 


B2'/(B + V + R) 
BT/(B + V + R) 
BO'/B 


+0.94 ±0.11 


+0.820 ±0.032 


+0.924 ± 0.022 


+0.15 ±0.05 


+0.099 ±0.025 


+0.150 ±0.029 


+0.65 ±0.13 


+0.449 ±0.068 


+0.705 ±0.063 



Table 5 
BSS Candidates 
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